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Abstract 

We are concerned with multidimensional stochastic balance laws driven by Levy processes. Using bounded 
variation (BV) estimates for vanishing viscosity approximations, we derive an explicit continuous depen¬ 
dence estimate on the nonlinearities of the entropy solutions under the assumption that Levy noise only 
depends on the solution. This result is used to show the error estimate for the stochastic vanishing viscosity 
method. In addition, we establish fractional BV estimate for vanishing viscosity approximations in case the 
noise coefficient depends on both the solution and spatial variable. 
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1. Introduction 

The last couple of decades have witnessed remarkable advances in the studies of partial differential 
equations with noise/randomness. A vast literature is now available on the subject of stochastic partial 
differential equations (SPDEs) and the particular frontier involving hyperbolic conservation laws with noise 
has had its fair share of attention as well. However, this is still very much a developing story and there 
still a number of issues waiting to be explored. In this paper, we aim at deriving continuous dependence 
estimates based on nonlinearities for stochastic conservation laws driven by multiplicative Levy noise. A 
formal description of our problem requires a filtered probability space (Cl, P, T, { r F,},> o) and we are interested 
in an L p (S. d )- valued predictable process u(t, •) which satisfies the Cauchy problem 

[du(t, x) + div v E(«(f, x))dt — [ n ri(u(t,x)\z) N(dz,dt), rellj, 

< ■'m > u (1.1) 

(u( 0, x) — uq(x), x € M. d , 

where = (R rf x (0, T)) with T > 0 fixed. The initial condition uq(x') is a given function on M. d , and 
F : R R' 1 ' is given (sufficiently smooth) vector valued flux function (see Section 0 for the complete list of 
assumptions). The right hand side of ( 11. lb represents the noise term and it is composed of a compensated 
Poisson random measure N(dz, dt) = N(dz, dt)—v(dz) dt, where A is a Poisson random measure on Rx(0, oo) 
with intensity measure v(dz), and the jump amplitude (integrand) rj(u,z ) is real valued function signifying 
the multiplicative nature of the noise. 

Hyperbolic conservation laws are used to describe a large number of physical phenomenon from areas 
such as physics, economics, biology etc. The inherent uncertainty in such phenomenon prompts one to 
account for the same and consider random perturbation of conservation laws. As an important first step 
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into the subject, a significant body of literature has grown around conservation laws that are perturbed by 
Brownian white noise. However, due to the complex nature of the uncertainties, it is only natural to look 
beyond Brownian white noise settings and consider problems with more general type of noise. We do that 
in this paper in the problem (11.11 ) by introducing Poisson noise in the right hand side. It is also mentioned 
the result of this paper could be extended to the general Levy noise case. 

In the case 77 = 0, the equation (11.11 ) becomes a standard conservation laws in R rf . For the deterministic 
conservation laws, well-posedness analysis has a very long tradition and it goes back to the 1950s. However, 
we will not be able to discuss the whole literature here, but only refer to the parts that are pertinent to the 
current paper. The question of existence and uniqueness of solutions of conservation laws was first settled in 
the pioneer papers of Kruzkov | 15] and Vol’pert 0 . For a completely satisfactory well-posedness theory 
of conservation laws, we refer to the monograph of Dafermos 10]. See also 0 and references therein. 


1.1. Stochastic balance laws driven by Brownian white noise 

As has been mentioned, evolutionary SPDEs with Levy noise has been the topic of interest of many 
authors lately, and new results are emerging faster than ever before. However, the study of stochastic balance 
laws driven by noise has so far been limited to equations that are driven by Brownian white noise and a 
satisfactory well-posedness theory is available by now. 

Observe that when the noise is of additive nature, a change of variable reduces equation into a hyperbolic 
conservation law with random flux which could be analyzed with deterministic techniques. In fact, Kim 
114] extended Kruzkov’s entropy formulation to establish the well-posedness of one dimensional stochastic 
balance law. 

However, when the noise is of multiplicative nature, one could not apply a straightforward Kruzkov’s 
doubling method to get a L 1 -contraction principle as in 111511 . The main difficulty lies in doubling the time 
variable which gives rise to stochastic integrands that are anticipative and hence the stochastic integrals in the 
sense of Ito-Levy would not make sense. Hence, it fails to capture a specific “noise-noise” interaction term 
relating two entropy solutions. This issue was first resolved by Feng & Nualart Ei with the introduction 
of additional condition, which captures the missing “noise-noise” interaction term, the so called strong 
stochastic entropy solution. They used L p framework to prove the multidimensional uniqueness result for 
strong stochastic entropy solution. However, existence was restricted to one space dimension since their 
proof of existence was based on a stochastic version of compensated compactness argument applied to 
vanishing viscosity approximation of the underlying problem. To overcome this problem, Debussche & 
Vovelle [9|] introduced kinetic formulation of such problems and as a result they were able to established the 
wellposedness of multidimensional stochastic balance law via kinetic approach. At around the same time, 
Chen & Karlsen [0] also established multidimensional wellposedness of strong entropy solution in L p D BV, 
via BV framework. Moreover, they were able to develop continuous dependence theory for multidimensional 
balance laws and, as a by product, they derived an explicit convergence rate of the approximate solutions to 
the underlying problem. We also mention that, using the concept of measure valued solutions and Kruzkov’s 
semi-entropy formulations, a result of existence and uniqueness of the entropy solution has been obtained 
by Bauzet. et. al. in JlJ]. 

In the article [11J], the authors used an entropy formulation which is strong in time but weak in space, 
which is in our view may give rise to problems where the solutions are not shown to have continuous 
sample paths. We refer to |3|], where a few technical questions are raised and remedial measures have been 
proposed. We also mention that Weinen et. al. IU8I1 published a very influential article describing the 
existence, uniqueness and weak convergence of invariant measures for one dimensional Burger’s equation 
with stochastic forcing which is periodic in x. 


1.2. Stochastic balance laws driven by Levy noise 

Despite relatively large body of research on stochastic partial differential equations that are driven by 
Levy noise, to the best of our knowledge, very little is available on the specific problem of conservation laws 
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with Levy noise. In fact, the first attempt were made to build a comprehensive theory on such problems 
in a very recent article by Biswas, et. al. J3l- For a detailed introduction to the SPDEs driven by Levy 
processes, we refer to the monograph by Peszat. et.al. IU6I1 and references therein. Roughly speaking, 
the theory developed in covers semi linear parabolic equations driven by Levy noise, which could be 
treated as stochastic evolution equations in some infinite dimensional Banach or Hilbert space, and typically 
the solutions of such equations enjoy regularizing properties. However, we can’t emulate those techniques 
on the specific problem of conservation laws driven by Levy noise due to the intrinsic discontinuous nature 
of the solution. 

In fact, independent of the smoothness of the initial data uq(x), due to the presence of nonlinear flux term 
in equation ( 11. 11 1. solutions to ( 11. Il l are not necessarily smooth and weak solutions must be sought. Before 
introducing the concept of weak solutions, we first assume that the filtered probability space (Q, P, 7L {J r t )t> o) 
satisfies the usual hypothesis, i.e., { r T t } t >o is a right-continuous filtration such that To contains all the P-null 
subsets of (Q, T). Moreover, by a predictable cr-field on [0, T] x O, denoted by Vj, we mean that the cr-field 
generated by the sets of the form: {0} x A and (s, t]x B for any A e Tq; B e T s , 0 < s,t <T. 

The notion of weak solution is defined as follows: 


Definition 1.1 (weak solution). An L 2 (lJ)-valued {Tt : 1 > 0}-predictable stochastic process u{t) - u(t, x ) 
is called a stochastic weak solution of (O if for all non-negative test functions if/ e C“([0, T ) x R rf ), 

I i/r( 0 , x)w( 0 , x) dx + I I x)u(t , x) + F(u(t, x)) • V v t/r(f, x)} dxdt 

Jr'' Jr j Jo ' 

+ I" r r T](u(t,x)\z)ip(t,x)dxN(dz,dt) = 0, P-a.s. (1.2) 

Ji =0 J|:|>0 Jr* 

However, it is well known that weak solutions may be discontinuous and they are not uniquely deter¬ 
mined by their initial data. Consequently, an entropy condition must be imposed to single out the physically 
correct solution. Since the notion of entropy solution is built around the so called entropy-entropy flux pairs, 
we begin with the definition of entropy-entropy flux pairs. 

Definition 1.2 (entropy-entropy fux pair). An ordered pair (ft, P) is called an entropy-entropy flux pair if 
P 6 C 2 (R) with p > 0, and y = (£i, fc, ••■•£/) : R i-» is a vector field satisfying 

£'(r) = p'{r)F'(r), for all r. 

Moreover, an entropy-entropy flux pair (/), P) is called convex if //'(•) > 0. 

With the help of a convex entropy-entropy flux pair (j 6, p), the notion of stochastic entropy solution is 
defined as follows: 


Definition 1.3 (stochastic entropy solution). An L 2 (R' / )-valued {'F, : t > 0)-predictable stochastic process 
u(t) = u(t, x) is called a stochastic entropy solution of (11. Il l provided 
(1) For each T > 0, p = 2,3,4, • • ■ , 

sup E[\\u(t, .)||£] < oo. 


0<l<T 


(2) For all test functions 0 < i/z e Cp ([0, oo) x Hr), and each convex entropy pair (J3 , P), 

I ip(0, x)P(u(0, x)) dx + I \dx)P(u(t, x)) + P(u(t, x)) ■ V x ifr(t, x)} dx dt 

Jr^ Jn r 1 

+ | I I ( P{u(r , x) + Tj(u(r, x); z)) - P(u(r , x))i//(r, x) dxN(dz, dr) 

Jr—0 J|’|>0 Jr^ V 2 

+ I I (p(u(r,x) + rj{u(r,x)\z)) - P(u(r,x)) - r\(u(r,x)\z)P' (u(r, x)))i//(r, x) v(dz) dr dx 
Jn T J| z |>o v 7 

>0 P-a.s. 
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Due to the nonlocal nature of the entropy ineaualities and the noise-noise interaction, the Definition 1 1.31 
alone does not seem to give the L 1 -contraction principle in the sense of average and hence the uniqueness 
is not immediate. In other words, classical “doubling of variable” technique in time variable does not work 
when one tries to compare directly two entropy solutions defined in the sense of Definion ll.3l To overcome 
this problem, the authors in |lj, y] used a more direct approach by comparing one entropy solution against 
the solution of the regularized problem and subsequently sending the regularized parameter to zero, relying 
on “weak compactness” of the regularized approximations. 

In order to successfully implement the direct approach, one needs to weaken the notion of stochastic 
entropy solution, and subsequently install the notion of so called generalized entropy solution (cf. Eli). 

Definition 1.4 (generalized entropy solution). An LrfW 1 x (0, l))-valued \'F, : t > 0}-predictable stochastic 
process v(f) = v(f, x, a) is called a generalized stochastic entropy solution of (11.11) provided 

(1) For each T > 0, p = 2, 3,4, ■ • • , 

sup E[||v(t, •,-)llp| < oo. 

0 <t<T 

(2) For all test functions 0 < i// e C' ,2 ([0, oo) x R 0 '), and each convex entropy pair (J3 , f), 

I" i/'IO, x)yS(v(0 ,x))dx+ I" I" (d t i//(t, x)j3(v(t, x, a)) + £(v(t, x, x)\dadxdt 

Jr^ Jn T Ja= o v 7 

+ | I II (/3(v(r, x, a ) + r/(v(r, x, a); z)) - J3(v(r, x, a))Wr, x) daN(dz,dr) dx 

Jr= 0 J|j|>0 JR^ Ja =0 V 7 

+ | | || (y8(v(r, x, a ) + rj(v(r, x, a)\ z)) - /3(v(r, x, a )) - rj(v(r, x, cr); z)y6'(v(r, x, a))j 

Jr =0 J|:|>0 Jr^ Ja=0 V 7 

x if/(r, x) da dx v(dz) dr 

>0 P-a.s. 

As we mentioned earlier, in a recent article J2l, the authors established well-posedness along with few 
a priori estimates for the viscous problem with Levy noise and proved the existence and uniqueness of gen¬ 
eralized entropy solution for multidimensional Cauchy problem (11.11 ) via Young measure approach. Finally, 
we mention that Dong and Xu 01 established the global well-posedness of strong, weak and mild solutions 
for one-dimensional viscous Burger’s equation driven by Poisson process with Dirichlet boundary condition 
via Galerkin method. Also, they proved the existence of invariant measure of the solution. 

1.3. Scope and outline of this paper 

The above discussions clearly highlights the lack of stability estimates for the entropy solutions of 
stochastic balance laws driven by Levy noise. In this paper, drawing preliminary motivation from (it, 
we intend to develop a continuous dependence theory for stochastic entropy solution which in turn can be 
used to derive an error estimate for the vanishing viscosity method. However, it seems difficult to develop 
such a theory without securing a BV estimate for stochastic entropy solution. As a result, we first address 
the question of existence, uniqueness of stochatic BV- entropy solution in UfW 1 ) n BViW 1 ) of the problem 
dm Making use of the crutial BV estimate, we provide a continuous depenece estimate and error estimate 
for the vanishing viscosity method provided initial data lies in uq 6 L p (R d ) n BV(R' / ). 

Finally, we turn our discussions to more general stochastic balance laws driven by Levy processes, 
namely when the function rj in the Levy noise term has explicit dependency on the spatial position x as well. 
In view of the discussions in 0, in this case we can’t expect BV estimates, but instead a fractional BV 
estimate is expected. However, that does not prevent us to provide an existence proof for more general class 
of equations in L p (R d ). 

The remaining part of this paper is organized as follows: we collect all the assumptions needed in the 
subsequent analysis, results for the regularized problem and finally state the main results in Section [2] In 
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Section [3] we prove uniform spatial BV estimate for the solution of vanishing viscosity approximation 
of (11. 1L and thereby establishing BV bounds for entropy solutions. Section 0] deals with the continuous 
dependence estimate, while Section 0 deals with the error estimate. Finally, in Section [6) we establish a 
fractional BV estimate for a larger class of stochastic balance laws. 


2. Preliminaries 


We mention that, throughout this paper we use C, K to denote a generic constants; the actual values of 
C, K may change from one line to the next during a calcuation. The Euclidean norm on any R^-type space 
is denoted by | • | and the norm in BV(S. d ) is denoted by | • 

Next, we collect all the basic assumptions on the data of the problem ( 11. 11 1 . 

(A.l) The initial function uq(x) is a np = i i 2 ,.T p (R' / )-valued (Fo-measurable random variable satisfying 

e[\\u q \\ p p + ||M 0 ||f + iMolffvpr)] < 00 for p - 1 , 2 ,.... 

(A. 2) For every k = 1,2..., d, the functions Fk(s) e C 2 (R), and Fk(s), F' k (s ) and F' k '(s) have at most polyno¬ 
mial growth in s. 

(A. 3) There exist positive constants 0 < A* < 1 and C > 0, such that for all u, v e R; z e R 

\B (h;z) - p(vz)\ < A*\u - v|(|z| A 1) 
and \rj(u\z)\ < C(1 + M)(|z[ A 1). 


(A.4) To prove existence and uniqueness of solutions, we assume that the Levy measure v(dz ) which has a 
possible singularity at z = 0, satisfies 



(1 A |z| 2 ) v(dz ) < +oo. 


Remark 2.1. Note that we need the assumption |(A.2)| as a result of the requirement that the entropy solutions 
satisfy L p bounds for all p > 2, which in turn forces us to choose initial data satisfying |(A.1)| However, 
it is possible to get entropy solution for initial data in L 2 (R' / ) n BV{ R' / ), provided the given flux function 
is globally Lipschitz. The assumption |(A.3)| is natural in the context of Levy noise with the exception of 
A* e (0,1), which is necessary for the uniqueness. Finally, the assumptions |( A. 1 )||( A.4)| collectively ensures 
existence and uniqueness of stochastic entropy solution, and the continuous dependence estimate as well. 

To this end, for any given fixed e > 0, we consider the viscous perturbation of (11. Il l 

du E (t, x) + di \ x F e {u e {t, x)) dt = I tj E (u e (t, x); z) N(dz, dt) + eA XK u e dt, t > 0, x 6 R d , 

4 l>o ( 2 . 1 ) 

m ( 0 , x ) = u e { 0 , x), x e R rf , 

where u e { 0 , x) is a smooth approximation of initial data uq{x) such that 

E[ f \u e (0,x)\ p dx] < £[ f \u 0 (x)\ p dx]. (2.2) 

M % M % 

Moreover, if initial data uq(x) 6 BV(R' / ), then 

J' |Vm £ (0, x)| dx J < J' |Vmq(x)| dx J. (2.3) 
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Furthermore, mainly to ease the presentation throught this paper, we assume that F e , r\ e are “sufficiently 
smooth” approximations of F and ij respectively. More specifically, we require that F, and ;/ f satisfy the 
same properties as F and q respectively (cf. |(A.2j| - |(A.3)| l and 

| F e (r) - F(r)\ < Ce(l + Ir^ 0 ), for some p 0 6 N, 

I Ve(u\z) - rj{u-,z)\ < Ce( 1 + |w|)(l a |z|). (2.4) 

Observe that, in view of J2, Subsection 3.2], these properties of I\ and q, are justified. 

For the deterministic counterpart of (12. It . proof of existence of global smooth solutions is classical by 
now. Same techniques could be used, mutatis mutandis, also for the stochastic scenario to establish the 
existence. More precisely, we have the following proposition from Jill . 

Proposition 2.1. Let the assumptions |(A. 1 )| |(A.2)| |(A.3)| and \(AA)\ hold and e > 0 be a given positive 
number. Then there exists a unique C 2 (M. d )-valued predictable process u e {t, •) which solves the initial value 
problem (12.1b . Moreover, 

(a) The solution u e (t, x) satisfies, almost surely, 

u e (t,x)~ I G(t, x - y)uo(y)dy - I I G(t - s, x - y)V ■ F e (u € (s,y)) dy ds 
JJs=o 

+ I I G(t— s,x-y)T/(u e (s,y);z)dyN(dz,ds), 

Js =0 J|j|>0 Jm$ 

where G(t, x) is the heat kernel associated with the operator eA xx i.e., 

1 

G(t,x ) G e (t,x) =- -e 4 « , t > 0. 

(4net)2 


( b ) For positive integer p = 1,2, 3, • • • , and T > 0 

sup sup £ , [||M e (f,-)llpl < 00 . (2.5) 

e>0 0<f<r 


(c) For a function fi e C 2 (M) with /3,j3' ,fi" having at most polynomial growth, 


sup E 

e>0 


e I I fi”(u e (t,x))\\/ X u e (t,x)^dxdt\ 

Jt =0 Jr 1 ; 1 


< 00, 


p = 1 , 2 ..., t > 0 . 


Remark 2.2. In view of Proposition ^. H and assumption |(A.l)| it follows that, for each fixed e > 0, Vufit, x) 
is integrable. Moreover if j RI , |V 2 « f (0, x)\ dx | < + 00 , then V 2 « f (f, x) is also integrable for fixed e > 0 and 
any finite time T > 0 (cf. |2| Section 3]). 

Now we are in a position to state the main results of this article. 

Main Theorem (continuous dependence estimate). Let the assumptions |(A. 1 )[ |(A.2)| |(A.3)| and |(A.4)] hold 
for two sets of given data ( uq , F, 77 ) and (vo, G, cr). Let u(t, x) be any entropy solution of (11.1b with initial 
data uq(x) and v(s,y) be another entropy solution with initial data i'o(y) and satisfies 

dv(s,y) + div y G(v(i,y))Js = I o-(v(s,y)-,z) N(dz, ds). (2.6) 

J|z|>0 
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In addition, we assume thatF", F'-G' e L“and define ©( 77 , cr) : = sup P -—-——— 

ueR J|z |>0 1 + [m|“ 

there exists a constant Cj > 0 , independent of |moIbv(r-') an d I v oIbv(R' , )> such that for a.e. t > 0 , 


v(dz). Then 


J' \u(t, x) — v(f, x)|0(x)Jx] 


< Cj 


(1 + £’[[vo|bv(R‘')]) y/t£>(ri, ct)||<X')|Il”(R‘') + £[|volw(R‘')] ll'f T/ ' ~ G’ IU ? ll0(’)ll/,“>(M rf > 



i«oW - v 0 (x)i <p(x)dx\ + 


(2.7) 


where 0 < cf) e C^(M. d ) such that | V0(jc)| < Cc/>( x) and |A0(x)| < C<p(x) for some constant C > 0. Moreover, 
a special choice of <p(x) with the above properties 


c/>(x) = 


(i. 

j e -c(ki-«) 


when |x[ < R , 
when \x\ > R , 


leads to the following simplified result: For any R > 0, there exists a constant C R r > 0, independent of 

l«olBV(R‘ i ) and IvoIbVCR-*). SUch that f ° r a ' e ' 1 - 


U | u(t, x) - v(t, x)| dx\ 

W</i J 

(1 + £[ | Vo | if VCR")]) y/tZXr],cr) + t E[|vo| B y(R<()] W F ’ ~ G 'IU + E \ f l«oM - v 0 (x)| dx] 

Jr i 


( 2 . 8 ) 


Remark 2.3. The condition that F", F' - G' e L°° could be avoided if we assume that u, v e L“((0, T) x 
K. rf x O) for any time 7’ > 0. In this case, an appropriate version of the main theorem would be possible. 
Moreover, the quantity ©(//, cr) is well defined in view of |(A.3)| and |(A.4)| 

As a by product of the above theorem, we have the following corollary: 

Main Corollary (error estimate). Let the assumptions |(A.1)| |(A.2)| |(A.3)| |(A.4)| hold and let u(t, x) be 
any entropy solution of (11.1b with E[\u(t, •)Ibv(r<')] < £[|koIbv(R‘ , )]» f° r t > 0. In addition, we assume that 
F" e L°°. Then, there exists a constant Cj > 0, independent of |koIbv(R‘'), suc l! that for a.e. t > 0 

E\ I \u e (t,x) - u(t,x)\dx\ 

L JrL 

< C T {e^(\ +£ , [|moIw(R‘')])(1 + t) + E\ J" |n e (0, x) - u 0 (x)\ </x]}. 

Moreover, if we assume that the initial error E [ |« f (0, x) - mq(x)| dx] = 0(e ^-), then we get 


E 


X \u e (t, x) - Ll(t, x) dx 

i 


0(eb 


Here we used the notation O(e) to denote quantities that depend on e and are bounded above by Ce, where 
C is a constant independent of e. 


Remark 2.4. We mention that, just like the deterministic case II12I1 . we are able to show that the rate of 
convergence for vanishing viscosity solution is i. It is also worth mentioning that this rate is optimal. 


7 

























We finish this section by introducing a special class of entropy functions which will play a crucial role 
in the analysis. Let f : R —» R be a C°° function satisfying 


m = 0, P'(-r) = -p'(r), 0" > 0, 


and 


For any f > 0, define /L : 


p'(r) = 
: by 


- 1 , 

e [-1,1], 
+ 1 , 


when r < — 1 , 
when |r| < 1 , 
when r > 1 . 


Then 




|r| -M x f< 0 ( {r) < |r| and |/£'(r)| < 

where denotes the characteristic function of the set A, and 

Mi = sup ||r| -0{r)\, M 2 = sup |y8"(r)|. 

W <1 kl<l 

Finally, by simply dropping for/? = we define 

Ff(a, ft) = JT /?V - ft)F» </(ct), F*(a, ft) = (ff (a, ft), f*(a, ft),ft)), 

F k (a, ft) = sign(a - b)(F k (a) - F k (b)), F(a, ft) = (Fi(a, ft), F 2 (a, ft), F d (a, ft)). 


(2.9) 


3. A priori estimates 


In this section, we derive uniform spatial BV bound for the stochastic balance laws driven by Levy 
process given by ( 11. Il l under the assumptions |( A. 1)||(A.2)|[(A. 3)[ and |(A.4)| 

Theorem 3.1 (spatial bounded variation). Let the assumptions |(A.1)| |(A.2)| |(A.3)| and |(A.4)| hold. 
Furthermore, let u e (t , x) be a solution to the initial value problem dp - Then, for any time t > 0 

|vw e (f, x)| r/xj < £■[ |VM £ (0,x)|dx] < e [ J' |vko(x)| nfx]. 

Proof. Since u e {t,x) is a smooth solution of the initial value problem (12. It . by differentiating (12.11) with 
respect to x,, we find that d Xi u e (t, x), 1 < i < d satisfies the stochastic partial differential equation given by 


d(d Xj u e (t, x)) + di y x (F' e (u e (t, x))d Xi u e (t, x)) dt 


I ij' £ (u e (t, x); z)d Xi u € (t, x)N(dz, dt) 

J[zl>o 


+ eA xx (d x u f {t, x)) dt. 


To proceed further, we apply Ito-Levy formula to ff:(d Xi u r (t, x)) to obtain 
d{fi c (d x u f (t, x))) + di v x (F’ e (uft, x))d x .u e (t, x))0^(d x .u e (t, x)) dt 
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lzl>o Je=o 


i)' £ {.u e {t, x); z)d Xi u e (t, x)P'^{d Xi u e (t, x) + 6 Tj' e (u e (t, x ); z)d Xi u e {t, x)) d6N(dz, dt ) 


(1 - 6){rj' e {u e \z)d Xi u e ) 2 Pf(d Xi u e {t, x) + 8rj' e (u e (t, x); z)d Xj u e (t, x)) d6 v(dz) dt 

|z|>o Je =o 

+ eA xx (d Xi u € (t, x))/^(d Xi w e (f, x)) dt. 

Since /3^ is convex, we conclude that 

eA xx (d x .u e (t, x))^(d Xi u e {t, x)) = e(Afe(d Xi u E (t, x)) - ^(d Xi u e (t, x))\Vd Xi u € (t, x)| 2 ) 

< eAfe(d Xi u f (t, x)), 

and for the martingale term, we have 

-l 

rj' E (u e (s, x); z)d Xi u E (s, x)[5'Ad Xi u e (s, x) + 8rj' e (u e (s,x)\z)d Xi u e (s, x))d8N(dz, ds ) 


(3.1) 


0 J|:|>0Je=0 


= 0. 


By Remark |2.2| we see that for each fixed e > 0 and 1 < i < d, S7d Xi u e (t,x) is integrable. Let 0 < ip(x) e 
C“(R' / ). Multiply (13.1b by if/ and then integrate respect to x to have 

J' /3^(d Xj u e (t, x))iA(x) dx j - E [ f / /3{(d Xi u e ( 0, x))iA(x) dx] 

(1 - 6){3”(d Xi u e (s,x) + 6ti' e (u e (s,x);z)d Xi u e (s,x)] 

x ( 77 ' (u € (s, x); z)d Xi u e (s, x))“i^(x) d8 v(dz) ds </x] 
div x (T E (u e (s, x))d x u e {s, x))/3^(d Xi u E (s, x))i//(x) ds r/xj 


< E 


i Js =0 J|:|>0 JS =0 


j=0 


+ eE 


fe(d Xj u e (s,x))Ai//(x) ds dx]. 


(3.2) 


To proceed further, observe that 

div x (F' e (u € (s, x))d Xl u e (s, x))^(d Xi u e (s, x))i//(x) = div x [ F'(w f (.v, x))d Xi u e {s, x)j3'^(d Xi u e (s, x))^(x)] 

- d Xi u e (s, x) F' E (u e (s, x))(p”(d Xi u E (s, x)) i//(x) Vd Xj u € (s, x) + f^(d x .u E (s, x)) Vi^(x)). 
Therefore, we obtain from (13.2b 

e[ f M d xM t,xmx)dx]^E[ f fe(d Xi u e (0,x))ift(x)dx] 

(1 - ff)/3g(d Xi u e (s,x) + 8t]' e (u e (s, x); z)d Xi u e (s, x)) 

x ( 77 ' (u e (s, x); z)d Xi u e (s, x))“( //(x) d8 v(dz) ds ofx] 

i I I d x .u e (s, x)ift(x)/3'J(d x u E (s, x))Vd x .u e (s, x) • F' e (u e (s, x))dsdx\ 

L Jr? Js =0 J 

i.C.f d Xi u e (s, x)/3^(d Xi u e (s, x)) Vi^(x) • F' e (u e (s, x)) ds t/xJ 

s I I P((d Xj u e (s, x))A(f/(x) dsdx\ 

L Jr? Js=0 J 


+ £ 


R? Js =0 J|z|>0 Js=o 

+ £[ 

+ £[ 

+ eE\ 
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(3.3) 


■ =e[ f f>d d x, u e(0, x))t]/(x) dx\ + £, (e, £) + £ 2 (e, £) + £ 3 (e, £) + S 4 (e, <f). 

To estimate &\(e ,£), we proceed as follows. Note that we can rewrite S\(e ,£) as 

£1 (e,£) = E\ f f f f (1 -0)h 2 p';(a + eh)il/(x)d0v(dz)dsdx\ 

1 Jr^ Js =0 J |:|>0 Js =0 J 

where a = d x .u e (s, x ) and h - q' f {u f (s, x); z)d Xj u e {s, x). In view of the assumption |(A.3)| it is easy to see that 


h 2 J3”(a + Qh) < |d. ri w e (s,x)| 2 (l A |z| 2 )/?"(a +0h). 


(3.4) 


Next we move on to find a suitable upper bound on a 2 /!'.'{a + 9 h). Since fi" is an even function, without loss 
of generality we may assume that a > 0. Then by our assumption |(A.3)| 

d Xi u e {t,x) + 0T]' e (u e (t, x); z)d x .u e (t, x) > (1 — A*)d Xj u e (t,x), 

for 6 e [0,1]. In other words 


0<a<(l-A*y\a + 6h). 


(3.5) 


Combining (13.4b and (13.5b yields 

h 2 %(a + dh) < (1 A |z| 2 )(l - AT 2 (a + 0h) 1 /3"(a + Gh)<C{ 1 A |z| 2 )f. 

Since by assumption |(A.4)[ (1 A |z| 2 ) v(dz ) < +°°, we infer that 

|£i(e,^)| < Cf£ m\^m and hence £i(e,£) i-> 0, as £ J, 0. (3.6) 

Next, we move on to estimate fiofe, F). In fact, we have 

\S2(e,&\<E ILL \d Xj u e {s, x)\>f/(x)fi"{d x M e {s, x))| \7d Xt u € (s, x)||£'(M e (s, x))| ds c/xj 

First observe that, in view of (12.9b , we obtain 

Mi 

| d Xj u e (s,x)\p'^(d Xi u e (s,x)) < \d Xj u e (s,x)\— X[-^4](d Xi u e (s,x)) h-> 0 , almost surely as F | 0, 
and moreover we see that 


\d Xi u e (s, x)| p'^(d Xi u e (s, x))tA(x)|v<9 Vi w e (s, x)| \F' e {u e (s, x))| 

< C||iA(0llt“(|Vi9x,Me(i,x)| 2 + |(M e (i,x))| 2/ ’°), for some po e N. 

In view of Remark 12.21 and Proposition 12.11 the right-hand side is integrable and independent of F > 0. 
Therefore, one can apply dominated convergence theorem to conclude that 

£ 2 (e^)w0,asn0. (3.7) 

Next, we consider the term £ 3 (e, £). With the help of uniform estimates (12.51) . we conclude 

|£ 3 (e,£)l < e\ f f \d Xi u e (s,x)\\Vi//(x)\\F' e (u e (s,x))\dsdx 1 
L Jr? Js= o J 
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(3.8) 


< I|V^-)IIz,”(]r.‘')£[ \d Xi u e (s,x)\\u e (s,x)\ P dsdx] 

jRi Js= 0 

< l|ViA(-)llL”(R‘ , ) £ ’f f f (\d Xi u e (s,x)\ 2 + \u e (s,x)\" P )dsdx] 

Jr? Js =o 

< C(e) T ||V^(-)|lz,»(Rrf), 


where we have used that for fixed e > 0 , d x u f {s, x ) is integrable. 

Finally we move on to estimate the term £ 4 ( 6 , £). It is easy to see that 

|£4(e,f)| < ||A^-)llL“(R‘') e£ f f f \d Xi u € (s,x)\dsdx] < 7’C(e)||A^(-)|lt«.(K") 

Jr? Js =0 

Taking advantage of (12.9b in (13.3b helps us to conclude 



x)|iA(x) Jx] < E[ f \d Xi u e (0, x)|i/r(x) dx] + ||tA(-)llz. 1 ( R rf ) 

yJ 

+ £i(e, 0 + £ 2 (e, £) + S 3 (e, 0 + £4(6,0. 


(3.9) 


(3.10) 


In what follows, we combine all the above estimates (13.6b . (13.7b . (13.8b . and (13.9b and then send £ i-> 0 in 
(13.10b to obtain 


i [ J' | d Xj u £ (t, x)\if/(x) dx] 

< E[ f \d Xi u e (0,x)\i//(x)dx] + C(e)(||A(KOIlL“ ( R«) + llV^OIb-oa*)) T (3.11) 


To this end, we define 0 < i/ov(x) e C^(R rf ) such that 


i// N {x) = 


when |x| < N 
when |x| > N + l. 


Note that since (13.1 lb holds for f{x ) = iI/n(x), we choose i//(x) = i/sn(x) in (13.1 lb . and then sending N —* 00 
to obtain 

£[ J \d Xi u e (t, x)| Jx] < Ey I \d X[ u e (0, x)| dx], 

which completes the proof. □ 


An important and immediate corollary of the uniform spatial BV estimate is the existence of B V bounds 
for the entropy solution of (11.1b . We have following theorem. 


Theorem 3.2 (BV entropy solution). Suppose that the assumptions |(A.2)| |(A.3)| and (AA)\ hold. Then 
there exists an unique entropy solution of (P > with initial data satisfying assumption |( A. 1 )\such that 


E[ I u(t, -)Ibv(r<J) ] < E[ \u 0 \ B vm ]’ f or an y f (3-12) 

Proof. We take advantage of the well-posedness results from and claim that the sequence {lif t, •)) con¬ 
verges, in the sense of Young measures, to the unique Z/(M d )-valued entropy solution u{t, ■). In view of the 
uniform BV estimate in Theorem l3.ll by passing to the limit, we conclude ( 13.12b . In other words, the unique 
//-valued entropy solution has bounded variation if the initial condition is BV. □ 
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4. Proof of The Main Theorem 

It is worth mentioning that, the average L 1 -contraction principle [see, for example, iH] gives the con¬ 
tinuous dependence on the initial data in stochastic balance laws of the type C£33 >- However, we intend to 
establish continuous dependence also on the nonlinearities, i.e., on the flux function and the noise coefficient. 
To achieve that, we need to consider the following regularized problem: 


idvfs,y) + diVyG e (v e (i,y ))ds = cr e (v e (s,y)\z)N(dz, ds ) + eA yy v e {s,y)ds, (s,y) e n r , 
{v e ( 0 ,y) = v*(y), yeW 1 ; 


where (v[L cr e , G,) are regularized version of (vo, cr, G) satisfying the conditions in ( 12.4b . In view of Theorem 
13.21 we conclude that v e (s,y) converges, as Young measures, to the unique BV-entropy solution i’(,v, v) of 
(12.6b with initial data vo(y). Let u(t, •) be the unique BV-entropy solution of (11. Il l with initial data uq(x). 
Moreover, we assume that the assumptions |(A.l)||(A.2)||(A.3)| and |( A.4)] hold for both sets of given functions 
(Vo, G, cr) and {uq, F, rj). 


We estimate the L 1 -difference between two entropy solutions u and v. The theorem will be proved by 
using the “doubling of variables” technique. However, we can’t directly compare two entropy solutions u 
and v, but instead we first compare the entropy solution u(t, x ) with the solution of the viscous approximation 
( |4~TT >, i.e., v f (,v, v). This approach is somewhat different from the deterministic approach, where one can 
directly compare two entropy solutions. For deterministic continuous dependence theory consult flSiQ 
and references therein. 

To begin with, let p and g be the standard mollifiers on M and respectively such that supp (p) c [-1,0) 
and supp(p) = B i(0). For 6 > 0 and > 0, let ps 0 (r) = —p(^) and g$(x) = jjg( |). For a nonnegative test 


function if/ e C*' 2 ([0, oo) x K (/ ) with x)| < C if(t,x), \Aif/(t,x)\ < Ci//(t,x) and two positive constants 

6, do, define 


x, s,y) = p So (t - s)g 6 (x - y)i//(s,y). 


(4.2) 


Observe that p$ 0 (t — s) £ 0 only if s - do < t < s, and therefore <t>s,g 0 (t, x; s,y) — 0 outside s - do < t < s. 

Furthermore, let g be the standard symmetric nonnegative mollifier on R with support in [-1,1] and 
gfr) = 7 ( 7 ( 7 ) for / > 0. We now write the entropy inequality for u(t,x), based on the entropy pair (J3(- - 
k), F^f, k)), and then multiply by gi{v e {s,y) - k), integrate with respect to s,y, k and take the expectation. 
The result is 


0 <e\ (ff fi(u(0,x) - k)c/)g'g 0 (0,x,s,y)gi(v e (s,y) - k)dkdxdyds\ 

1 Jn T Jr? Jr 4 ’ 1 

: I I /3(u(t,x) - k)d,(pgg 0 (t, x, s,y)gi(v e (s,y) - k)dkdxdtdyds\ 

L Jn r Jn T Jrj ’ J 

((( (j3(u(t, x) + r/(u(t, x); z) - k) - J3(u(t, x) - k)) 

i\.j J Jn r vj|2|>o 

x <, x, s, y) gi(v f (s, y) - k) N(dz, dt) dx dk dy c/s] 

■£[ r r r r r m xy, z) ~ k) ~piuf, x) - ^) 

L Jn r Jt= 0 J|zi>o Jr? Jri v 


Jn T * 

+ E\ 

+ E 


- rj(u(t, x); z)/3’(u(t, x) - k 


, x; s,y) 


x gi(v e (s, y) - k) dk dx v(dz) dt dy c/sj 
E\ f f f F fi (u(t, x), k) ■ V x Qs(x - y) f(s,y)p 6o (t - s) 

v/flr v/Ilr '-'Mt 
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x Q(v e (s,y) — k) dkdxdt dy afsj 


I\ + 12 + h + U + h- 


(4.3) 


We now apply the Ito-Levy formula to ( 14.1b and multiply with test function (f>g 0 ,s and $i(u(t, x) - k ) and 
integrate . The result is 

0 <Ilf j8(v e (0,y) - k)<pg'g 0 (t,x,0,y)gi(u(t,x) - k)dkdxdydt\ 

1 Jn T JrJ Jr* ' " J 

+ E II I I /3(v e (s,y)-k)d s (f>gj 0 (t,x,s,y)si(u(t,x)-k)dkdydsdxdt\ 

1 J n r Jn r Jr* ’ J 

+E \ f f f f (j 8 ( v e(s,y) + o- e (v e (s,y);z)-k)-jB(v e (s,y)-k)) 

*^117' Jn T *-/|z|>o 

x <f>gj 0 (t, x, s,y)gi(u(t, x) — k) dkN(dz, ds) dy dxdt] 

+£ t f f f f f (yS(v f (.v,y) + cr f (v f (.v,y);z) - fc) - [i{v e (s,y) - k) 

Jn T Js=o J| z |>o Jr^ Jr* 

- cr e (v e (s,y); z)P'(v e (s,y ) - k))<p Sy6o (t, x; s,y ) 

x Q(u(t, x) - k ) dk dy v(dz) ds dx dt j 

G^(v e (s,y), k) ■ V y Qg{x - y)i]/{s,y)pg a (t - s) qi{u{t , x) - £) dkdxdtdyds J 

Gf (v e (s, y), £) ■ V y i//(s , y)ga(x - y)pg 0 (t - 5 ) Si(u(t, x) - k) dk dx dt dy c/sj 

- e£j | | f y 6 '(v e (s,y) - &)V y v e (s,y) ■ V y (f>s,g 0 (t,x, s,y) $i(u(t, x) - k)dkdy ds afxoffj, (4.4) 

*-'117' 

where (a, b) = _£* fi'(r - b)G' 6 (r) dr. It follows by direct computations that there is /? € N such that 

|Gf(a, b) - G^(a, fe)| < Ce(l + |a[ 2p + |fe| 2p ). 

In view of the uniform moment estimates, it follows from (14.4b that 


+ E 

+ E 


p lip H&./, 


J kJ Up ]RL£ 


0 <E 


n r JRi Jr* 


y6(v e (0, y) - k)<f>g'g 0 (t, x, 0, y)$i(u(t, x) - k) dk dx dy Jf| 


+ E 


Jiip Jn T 


+E 


117- *^117- *j|z|>0 *JRj. 


j 8 (v e (s, y) - k)d s c/>g y g 0 (t, x, s, y)pi(u(t, x) - k) dk dy ds dxdt j 
f (jS(v e (s,;y) + cr e (v e (s,y);z) - k) -p(v e (s,y) - k)) 

JR k 

x (ps t 5 0 (t, x , s , y)$i(u(t, x) — k) dkN(dz , ds) dy dxdt\ 

+£ t f f f f f (j 3 ( v e( s >y) + °'e(v e (s,y);z)-k)-p(v e (s,y)-k) 

Jn T Js=o J|j|>o Jr^ Jr* 

- cr e (v e (s,y);z)P'(v e (s,y) - kj)(f> SySo (t, x; s,y ) 

x $i(u(t, x) — k) dk dy v(dz ) ds dx Jf] 

G^(v e (s, y), k ) • VyQgix - y)ip(s, y)pg 0 (t - s)si(u(t, x) - k) dk dx dt dy J.vj 


+ E 


p lip kj 
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+ E 


-eE 


m T Jn T 


G^(v e (s, y), k) ■ V y tJ/(s, - y)ps 0 (t - s)$i(u(t, x ) - k) dk dx dt dy Js] 

P'(v e (s,y) - k)S7 y v e (s,y ) • ^y<t>s,s n Q(u(t, x) - k ) dkdydsdxdt^ + C(J3, iff)— 


/[ + J 2 + J 3 + J 4 + J 5 + A + Ji + C(J3,i//)~, 

o 


(4.5) 


where CY/i, 1 //) is a constant depending only on the quantities in the parentheses. Our aim is to add (14.3b and 
(TO . and pass to the limits with respect to the variousparameters involved. We do this by claiming a series 
of lemma’s and proofs of these lemmas follow from 12 |] modulo cosmetic changes. 

To begin with, note that particular choice of test function (14.21) implies that Jj = 0. 

Lemma 4.1. It holds that 


h + 


J\ —> E\ f f f f)(u(0,x) - k)ij/(0,y)g s (x - y)Q(v f (0,y) - k) dkdxdy\ 
^o-* 0 JjRrf jRi J Ki 

7 —* e\ r f J3(u(0,x)-v e (0,y))i/r(0,y)g s (x-y)dxdyl 
L Jr^ Jr; j 

We now turn our attention to ( I 2 + Ji)- Since /j, 57 are even functions, we see that 


h + Ji 


= 4 j f f 


P(v e (s,y) - k)d s i//(s,y)p So (t - s)g s (x - y) 

k 

x gi(u{t,x) — k) dk dy ds dxdtf 

Lemma 4.2. It holds that 

h + Ji —> E\ I I I / 3 (v € (s,y) - k)d s i/j(s,y)gs(x - y)$i(u(s,x) - k) dk dy dx ds\ 

L J n ^ J Kj ' J 

-j-^e[J J' J 3 (v e (s,y) - u(s,x))d s ft(s,y)g s (x-y)dydxds]. 

Next, we consider the term Is + Js and regarding these terms we have the following lemma. 
Lemma 4.3. The following hold: 

! II F^(u(s, x), v e (s,y)) ■ V x gg(x - y) if/(s,y) dxdy ds\ 

L Js =0 Jr; Jr; j 


lim lim Is - E\ 

/-»0 fi o -»0 


(4.6) 


and 


lim lim Js = e\ f f f G^(v e (s,y),u(s,x)) ■ V v g#(x - y) i/r(s,y) dxdyds] (4.7) 

1^0 60^0 L J s =oJr;Jr; ' 1 


Lemma 4.4. It holds that 

J(, —> E l| | | G^{y e {s,y),k) ■ V y f{s,y) gs{x - y)gi{u{s,x) - k) dkdxdy ds\ 

*0 L J n ^ J R ^ .... . j 

—> £ | | G /3 (v € (s,y),u(s,x)) ■ V v iff(s,y)gs(x - y) dxdy ds\. 

(-*o L Jn T Jr; j 
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Next, we consider the term J 7 . Thanks to the uniform spatial BV estimate for vanishing viscosity solution 
(cf. Theorem 13. 1 1 . we conclude that 

\Jl\ < eP'lloo^f f fjv,v^)||v^.^- y) |^]| 
v/ n r 

< e\\P'W°c,E\ f f f \VyV e (t,y)\\Vy[i//(t,y)g 6 (x-y)]\dxdtdy] 

J\ y \<K Jt =o Jr? 

^ C^^IIvoIbvcr'')] (4.8) 

Lemma 4.5. It holds that 

lim lim J 4 = e\ f f f f (1 - A)j8"(v e (s,y) - u(s,x) + Acr e (v e (s,y);z)) 

i -* 0 Jn r Jr? J| Z |>o Ja =o 

x \o- e (v € (s,yy,z)\ 2 ilj(s,y)Qs(x - y) dAv(dz) dxdy ds], (4.9) 
lim lim U = C f f f (1 - A)/3"(u(s, x) - v' e (s,y) + Ai)(u(s, x); z)) 

x \rj(u(s, xy,z)\ 2 ^(s,y)gs(x - y) dAv(dz) dxdy ds\. (4.10) 

Finally, we consider the stochastic term / 3 + J 3 ; 

Lemma 4.6. It holds that J 3 = 0 and 

lim lim / 3 = e\ I I I (/3(«(r, x) + r](u(r, x); z) - v e (r, y) - cr e (y £ (r, y); z)) 
i-> 0 L Jn,- Jr? J| z |>o v 

~/3(u(r,x) - v E (r,y) - cr e (v e (r,y);z)) +/3(u(r,x ) - v e (r,y)) 

-P(u(r, x) + rj(u(r, x); z) - v e (r,y))jtfi(r,y) g e (x - y) v(dz) dxdydr^. 

To proceed further, we combine Lemma l4~6l and Lemma H31 and conclude that 


lim lim f (/ 3 + J 3 ) + (I 4 + J 4 )) 

/^o <5 0 -^o v ’ 

=e\ f f ( f {/3(u(t,x)-v £ (t,y) + 7j(u(t,x);z)-o- £ (v £ (t,y);z)) 

Jn T Jr? J|z|>o 

- j8(u(t, x) - v e (f, y)) - (77 (u(t, x); z) - cr e (v e (t, y); z)) 
x P\u{t, x) - v e (f,y))} v(dz)^i//(t,y)gs(x - y) dxdy c/fj 

=£[ f f f f f /3"( u ( r >x) — v e (r,y) + p(r/(u(r,x);z) - cr e (v e (r,y);z))) 

J r =0 J|zi>o Jr? Jr? J p =o 

x (1 - p)\rj{u(r, x); z) - cr e (v e (r, y); z)|>(r, y) 

x p^(x - y) dp dx dy v(dz) dr J (4.11) 

We are now in a position to add (14.3b and (14.5b and pass to the limits lim lim. In what follows, invoking 

/—>0 dotO 

Lemma |4J~1 Lemma l4~2l Lemma l4~3l and Lemma l4~4l and the expressions (14.8b and (14.1 lb . we arrive at 



u( 0, x) - v e (0, y))if/( 0, y)gs(x - y) dx c/yj 
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+ E 


-[/ / 

Jr =0 J|z|>0 


[ P(v e (s,y) - u(s,x))d s i//(s,y)g s (x-y) dydxds] 

Jn T Jr? 

- e\ I I V y ■ {G^(v e (s,y),u(s,x)) — F^(u(s,x),v e (s,y))}i//(s,y)gg(x—y)dydxds] 

1 Jn T J r? 1 

?[J J'^E /3 (u(s,x),v e (s,y)) ■ V y ifr(s,y)g d (x-y)dydxds] + c(£[|v 0 | BV(R <i)] + l)^ 
P"(u(r,x) - v e (r, y) + p(rj(u(r, x); z) - cr e (v e (r,y);z ))) 


i jp =o 


x (1 - P)\tj(u(r,x)\z) - o- e (v e (r,y);z)|Xr,y) 

x gs(x - y) dp dx dy v(dz ) Jr] 

!J{\ + + ^3 + + X + C^£’[|voIbv(R'')] + l)~- (4.12) 


Again, our aim is to estimate all the above terms suitably. First observe that, since p^(r) < |r|, we obtain 

1X1 < £■[ f f |v e (0,y) - u{0,x)\\l/{0,y)g s (x-y)dxdy\. (4.13) 

My JM.^ 

Next, by our choice of p = fy, we have 

| |;(F^(«, v) - F*(y, u))\ = | - F\v)^(y - u) - F'(v)%( 0) + J j8"(s - v)F'(j) ds| 

- |(F'(v) - F'(u))%(u - v) - f - v)F"(j) dsl 

J S-U 

= I f (/?'(« - v) - /?'(s - v))F"(i) *1 < M 2 f ||F"|U 

Jit 


Also from the definition of F& and G&, it is evident that 

||;(^(v, u) - GP(y, m))| < |F'(v) - G'(v)| 
Therefore, by ( 14. 14b and ( 14. 1 5b . we obtain 

|^-(f^(m,v)-G^(v, h))| < M 2 %\\F"\U + | F\v) - G'(v)| 
Keeping in mind the estimate (14. 16I ). we proceed further by rewriting the term X as 


(4.14) 


(4.15) 


(4.16) 


X =£ 


[ f f V v v e (s,y)-5,,(F' /J (M,v)-G /J (v,M))| ^(s,y)^(.r-y)£/y<lrr/.s| 

Jn r Jr-* K«.v)=(«(*,*),v e (ij)) J 


Thanks to the uniform spatial BV estimate for vanishing viscosity solution (cf. Theorem l3. II) . we conclude 
that 

\W 3 \<(M 2 £\\F"\\ 00 + \\F'-G'W^E i JILL | V y v c (s, yM(s, y)g s (x - y) dx dyds\ 


^ £’[|voIbv(R'')]( 4^2 £ l|T’ ,, ||oo + ||F' - G'Hoo) ( ll^(s, Olb-cn*)^- 

J s =0 


(4.17) 
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Next, we recall that the function i/c(f, x) satisfies \Vi//(t, x)| < C t/c(f, x) and \F^(a, b)\ < ||F'||oo|cj - b\ for any 
Therefore, we conclude 


\9U\ <C[|F'|| l ~f[ f f f |m(.s\ x) - v f (.s\ y)\if/{s, y)p,s(x - >■) dxdy ds\ 

^ = 0 My 

<C||F'|| l ~f[ f f f J3f(u(s,x)-v e (s,y))i/r(s,y)p#(x-y)dxdyds 1 
Jj=0 Jr? Jr? 

f i 

Js=0 


+ CM, ||F'| 


ll<A(s> •)[[l”(R' j ) Js- 


(4.18) 


Let us focus on the term JL. For this, let us define 

a := u(r, x) - v e (r,y), and b q(u(r, x); z) - <x e (v e (r,y); z). 
Then JL can be rewritten in the following simplified form 


#5 


= £ i LL 
[£X 


" Jp=0 


( 1 - p)b^ft" (a + p b)i//(r, y) p^(x - y) dp dx dy v(dz) dr j 


< C£ 


? j p = 0 


|/ 7 (i/(r, x); z) - cr(ic(r, x); z)| ft"{a + pb) 


x i //(r, y) p,;(x - y) dp dx dy v(dz ) c/rj 

' C£ t f f f f f \o-(u(r,x);z)-o-(v f (r,y);z)\~j3"(a+pb) 

J r = 0 J|zi>o Jr? Jr? Jp=o 

x i/c(r, y) £>,;(x - y) dp dx dy v(dz ) c/rj 


+ CE 


^ io X l>o X? 


5 ^p=o 


|cr(v e (r, y); z) - <x f (v f (r, y); z)|“j8"(a + p b) 


X i/c(r, y) g^(x - y) c/p c/x c/y v(c/z) c/rj 


:= ^5 +^5 


(4.19) 


To this end we recall that £>( 77 , cr) = sup 

t/eR J|z|>0 

|(A.3)| With this quantity at hand it is easy see that 


■f 

: J|~l= 


|/7(m, z) - ct(m, z)|* 

1 + \u\ 2 


v( c/z), which is well-defined in view of 


^5 < 


CTXq, cr) 


CD(rj, cr) 


U ( 1 + |«(r, x)| 2 )i/c(r, y)p rf (x - y) c/y c/x c/rl 

r=o Jr? Jr? 

(X II^(s>')[Ili + ll<A( r , Olloot/r) 


(4.20) 


Next, we move on to estimate the term ??[ 2 . Observe that 


|cr(«(r,x);z) - o-(v E (r,y)-z)\~P"(a + pb) < \u(r,x) - v e (r,y)| 2 (l A \z\ 2 )/3"(a + pb) 

= (1 A\z\ 2 )a 2 p"(a+pb). (4.21) 


Therefore, it is required to find a suitable upper bound on crfi"(a + pb). Since //" is non-negative and 
symmetric around zero, without loss of generality, we may assume that a > 0. Then, by our assumption 
|(A.3)[ we conclude that 

\q{u(r, x); z) - cr e (v e (r, y); z)| 
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< 177(w(r,x);z) - cr(u(r,x);z)\ + \cr(u(r, x); z) - cr(v e (r,y); z)| + |<x(v e (r,y);z) - cr e (v e (r,y); z)| 

< \r](u(r, x);z) - cr(u(r, x)\z)\ + A* a + Ce( 1 + |v e |), 

which implies that 

a + pb > -\r/(u(r, x);z) - cr(u(r, x);z)| - Ce( 1 + |v e |) + (1 - A*)a, 
for p 6 [0,1]. In other words 

0 < a < (1 - /l*) _1 {fl + pb + |t 7 (m(t\ x);z) - cr(u(r, x);z)| + Ce( 1 + |v f |)}. (4.22) 

Now, we shall make use of ( 14.221 ) in ( 14.211 1. to obtain 

|cr(M(r,x);z) - cr(v e (r,y);z) \ ( a+pb) 

< (1 - /fT 2 {(fl + pb) 1 + C \r](u(r, x)\z) - cr(u(r, x)', z) f + Ce 2 (l + |v e | 2 ) }(1 A \z\ 2 )^(a + pb) 

1 A |z| 2 ). 

This helps us to conclude 


< c(% + ~ ^“( r ’X);z) r + e2 1 1 + 


I^g| < CE f f f f U+ —-— --)(1 A \z\ 2 )il/(r,y)g s (x-y)dxdym(dz)dr 

Jr J| Z |>0 Jr? Jr? v £ ' 

rr r 

Jo Jr? Jr? 


£>(r1, cr ) 


(1 + \u(r, x)| )i]/{r,y)ps{x - y) dxdy dr 


<C(t + 


7") f Ilf 0, ■)llz.°°(R rf ) Js + C!D ^’ ( f ||f(,s, Olli 1 ds+ f llf(r, OlLJr). (4.23) 

S Js =0 S Jo Jo 


Next, we move on to estimate the term JT 2 In fact, it follows that 


E\ f f f f f \cr(v, ! (r,y);z)-cr 6 (v e (r,y);z)j 2 j3"(a+pb) 

Jr—O J|z|>o Jr? Jr? J p = 0 

X y) g$(x - y) dp dx dy v(dz) Jr) 


,C E [f f f f 

Jr= 0 J|:|>0 Jr? Jr; 


! (l + |vj 2 ) 


(1 A \z\ L )^(r,y)Q S (x-y)dxdyv{dz)dr\ 


A rT 


< C- 


f 1 


lifts, 0 IIl»(r-A ds. 


(4.24) 


We now make use of the estimates (14.20b . (14.23b and (14.24b . Then it is evident from ( 14.19b that 




CD(Tj,cr) 




II f(s, 01b ds + 


£ IlfOIL Jr) 


+ C(t + 


f Ilf ( s , 0 IIl” ( r") ds. (4.25) 

? Js=0 


Finally, we make use of the estimates (14. 13b . (14. 17b . (14.18b and (14.25b in (14.12b and pass to the limit as e —> 0 
(keeping 6 and j fixed) in the resulting expression to conclude that 


0 <E 


r? Jr? 


|v' 0 (y) - m( 0, x)|i/r(0, y) Q S (x - y) dxdy\ 


18 

































■ £’[kolBV(K' i )](^2^|[^’ ,, |[oo + \\F' - G'lloo) r l|t/ r (i, •)IIl“(R-') ds 

J 5=0 

■ C\\F'\\ l ~e\ f f f fe(u(s,x)-v(s,y))i//(s,y)Q S (x-y)dxdyds\ 

Js =o Jr? Jr? 

+ C(M\\\F'\\ij*. + V)% r •)IIl“(r-')^'S 

Js= o 


CD(rj, cr) 


(X iw "’- )iiiijs+ X ||iA(r, Olloo^r) 

I I /?(v(s,y) - m(s, x))d s ifr(s,y)gs(x - y) dy dxds . (4.26) 

L J n T Jr? j 


Now we can safely pass the limit as 5 —> 0 in (14.26b to obtain 
0 <ft[ J' |vq(x) - u( 0, x)|i/r(0, x) dx j 


■ £’[| m oIbv(R‘o]( 442£||.F’ ,, |U + ||F' - G'Hoo) f 110'C-s', •)llz.“(M rf ) ds 

J s= 0 

+ C||F'||£»£f f f J3f(v(s,x) — u(s,x))ifr(s,x)dxdsl 
L J s =o Jr? j 

+ C(Mi [|F'||£~ + 1)£ I" ||i/f(s, •)IIl“(r- j ) Js 

Js=0 

+ CTKi], cr) ^ ||^( iS) .)|| tl ds + ||^(r, -)l|oo dr) 

+ £j J' (3f(u(s,x) - v(s,x))d s il/(s,x)dxds]. 


(4.27) 


To proceed further, we make a special choice for the function x). To this end, for each h > 0 and fixed 
f > 0, we define 


= 


1, if s < t, 

1 — if t < s < t + h, 

0, if 5 > r + h. 


Furthermore, let <p 6 C^(K d ) be a cut-off function such that |V0(x)| < Ccp(x), |A0(x)| < Ccf>{x). Clearly, 
(14.271) holds with i//(s,x) = if/‘ h (s)cf>(x). Let T be the set all points t in [0, oo) such that t is right Lebesgue 
point of 

A(s) = _fi| J' j3f(v(s, x) - u(s, x)) cf>(x) dx J. 

Clearly, ^(complement of T) has zero Lebesgue measure. Fix t e T. Then from (14.271 ), keeping in mind 
that we used generic /3 for the function (3^, we obtain 

0 <£[ J' |t ; o(x) - u( 0, x)|i^(0, x) dx J 

+ £’[|' ; oIbi/(R‘')](442 £ ||F"|| 0O + || F’ - G'|U)|W0IIl-(*) f i/^(s) ds 

J s-0 

+ C\\F'\\ l ^e\ f f Pc{v{s, x) - u(s, x))if/‘ h (s)<p(x) dxds] 

Js=o Jr? 
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(4.28) 


CD(rj, cr) 


( f f (f>(x)<l/‘ h (s)dxds + f if/' h (s)\\(l)\\oodr) 
Jo Jr-' Jo 


— — f e\ f J3((u(s,x) - v(s,x))</>(x)dx]ds. 
h J s=t 1 J-Rd J 

Since t is a right Lebesgue point of A(s), letting /; —» 0 in (14. 28l > yields 


i[ J' j8f(u(t, x) - v(t, x))(f>(x ) dx j 

<£f f |v 0 (x) - u(0,x)\ip(x) dx] + C{M\ ||C'|| L «. + 1)£||<K-)IIz.”(r-') t 

Jr? 

+ £[|voI W (r-)](m 2 ^||F"|[oo + - G , || 0O )||^(-)|| i .( R , ) t 

+ CWF'Wu* f £[ I" /3g(v(s, x) - u(s, x))<p(x) dx\ ds 
Js=o L Jr? j 


CtD(j), cr) 


(||0IIl. + llflb-) 


for almost every t > 0. An weaker version of Grownwall’s inequality then yields 

E [ f R A (u(t ’ x) _ v(t ’ xmx) dx \ - eC, " F '"°° E [ £„ |v'o(x) - u( 0, jc)|0(x) dx j 

+ Ce C||p|l ” J j(M l HF'lli- + 1 )£||<K')IIl“(R‘' ) t 


+ £’[|voIbv(r-')](4J2 %\\F"\\oo + \\F' - G'||oo)||<^(-)[Il“(r-') t 
+ £^ (lwlil + lwll .)] 


(4.29) 


for almost every t > 0. Next, we recall that |r| < /f f (r) + M\ J, for any re R. Using this inequality, ( 14.29b 
reduces to 

E\ f f3 e (u(t, x) - v(t, x))</>(x) dx 1 < e c ' l|r|l ”£[ f |v 0 (x) - u( 0, x)\<p(x) dx 1 + M x ^||0(-)IIl 1 (r-') 

Jr? Jr? 

+ Ce c " F '"°° j(M! ||F'|b- + 1)£||0 (-)IIl"(rV 
+ -E[|v 0 Iw ( r-')](^ 2 f ||F"|L + II F' - G'|| 


CtD{T),(r) 


(ll0lli> + II0IIL-) 


(4.30) 


We now simply choose J = sJtTXq, cr) and conclude that for a.e t > 0 

E[£y{ t ,x)-v{t,x)\^x)dx\ <C T e[ £ d I u Q (x) - voMI <f>(x) dx\ + Sflvobv^)] II F' - G'|U t ||0(-)llz.“(R")] 

+ Ct [(1 + £’[|voby(R-')]) y/ED(rj, o")ll0(-)lli~(R rf ) + y/t!D(ri, o')II0(-)IIl 1 (R'')]> 

for some nonnegative constant Cr, independent of |moIbv(r-') and |v'olsv(R rf )- This completes the first part of 
the proof, and second part follows from this by exploiting the specific structure of the test function <!>(x). 













5. Proof of The Main Corollary 


It is already known that the vanishing viscosity solutions converge (in an appropriate sense) to the unique 
entropy solution of the stochastic conservation law. However, the nature of such convergence described by 
a rate of convergence is not available. As a by product of the Main Theorem, we explicitly obtain the rate of 
convergence of vanishing viscosity solutions to the unique BV-entropy solution of the underlying problem 

(TO . 

By similar arguments as in the proof of the Main Theorem (cf. Section^, we arrive at 


+ E 


0<£[ f [\u e {Q,y)-uo{x)\if/{Q,y)es{x-y)dxdy\ 

Ky JRJ 

i [l M ol,By(R‘o]Af 2 £ ll-F^IL ll<A(s, •)llt“(R'') ds + C— J' ||tK s > Ollz,”^) ds 

+ C\\F'\\ I ~E\ fe(u e (s,y) - u(s,x))i//(s,y)g s (x-y)dxdyds] 

Js=o Jr? Jr? 

+ C(M l IIF'Hi. + lJf J ■)IIl“(R‘') ds + C(l + £ , [|moIbv(r<')])- 

i I I J3/:(u £ (s,y)-u(s,x))dsi//(s,y)Qg(x — y)dydxds\. 

L Jn r Jr? j 


+ E 


(5.1) 


Let iff(s,y) = 4t' h {s)(p{y) where i/f h (s) and <fi(x) are described previously. Let T be the set all points t in 
[0, 00 ) such that t is right Lebesgue point of 

B(s) = E\ \ Pt;(u e (s,y) - u(s, x))<f>(y)g e (x - y) dxdy 1. 

Jr? Jr? 

Clearly, T c has zero Lebesgue measure. Fix t e T. Thus, from (15.11 . we have 

7 f E\ f [ /3 ( (u e (s,y)~ u(s,x))<f>(y)g e (x-y)dxdy]ds 
n J,=t L Jr? Jr? j 

r*t+h r* r* 

< C\\F'\\ l *, ( E \ \ \ <t>(y)fe(u e (s,y)-u(s,y))e s (x-y)tjf h (s)dxdy\ds 

Js =0 Jr? Jr? 

*LL |we( 0 ,y) - u 0 (x)\cf>(y) Q S (x - y)dx dy\ 

+ C£ , [|MolBV(R<')]Af 2 ^||T' ,, |UII<^( 0 llL“(R-') f ^/ h {s)ds + C — ||0(-)|[i~(R</) f t// h (s)ds 

Js=0 < Js=0 

+ Qr II<K p )IIl“(r-o t l / l(s)ds + C(l + £ , [|moIbv(R‘ j )])^- 


+ £ 


Taking limit as h —» 0, we have 
E 


\ Pc{u e (t, y) - u(t, x))<p(y)p s (x - y) dxdy] 

L Jr? Jr? j 

C\\F'\\v° f E\ f f </>(y)j3f(u e (s,y)-u(s,y))pa(x-y)dxdylds 
Js =0 Jr? Jr? 

M f (0,>’) - u 0 (x)\<f>(y)g s (x-y)dxdy\ + c(l + EOkoIbvtwo])^ 
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+ 


c(l + £[|^|fly(R<^])f ||0 (OIIl”(R<<) t + C — ||0(Ollz."(R'') t 


By an weaker version of Gronwall’s inequality, for a.e t > 0 

E\ \ P((ue(t, y ) - u(t, x))</>(y)Qs(x - y) dx dy 1 

L JR? JR? J 


< gC||F'|| t «o r 


\u e (0,y) - u Q (x)\cf>(y)g s (x-y)dxdy] + c(l + £[|m 0 |bv(r</)])^} 


■ Ce C + £[|MolBV(R‘ i )])^ IW0llz.“(R' i ) f + — H0(')llL”(R‘ i ) * 


Since |r| < M \<f + y0f(r), we have 
E 

„C\\F‘ 


ULJj MfCky) - u(t, x)\(Ky)Qs(x - y) dxdy ] 

w e (0,y) - M O (x)|0(y)£) I 5(x->')r/xdy] + c(l + £[|m 0 |bv(r<')])|} 


+ Ce c " F " i “'^l + £[|moIbv(R'')])^II^( - )IIl”(R‘') t + — II<^(-)IIl“(R'') f ] + C£ ||0(-)llz.”(R‘ i ). (5.2) 

First we send (f> to xk in (15.2b . and then choose £ = e. The resulting estimate gives 


f \u e {t,y) - u(t,x)\g s (x-y)dxdy\ 

L Jr;! Jr? 

'{e\ \u € (0,y)-u 0 (x)\g s (x-y)dxdy]+C(l + £ , [|m 0 Ibv(R‘')]) 

Jr? Jr? 


R? Jr? 

.C||F'|li' 


. (j e c H F Ik" '^1 + £ , [|MolBy(Rii)])e t + Ce. 


Notice that, 
E 


[ JJu e (t,y) -U(t,y) I dy] 

<£1 I I \u e (t, y) - u(t, x)\gs(x - y) dx dy + E\ I I \u(t, x) - u(t,y)\gs(x - y) dxdy\ 

L JR? JR? - J L J R j - - J 

<e\ \u e (t,y)-u(t,x)\g s (x-y)dxdy\+6E\\u 0 \ Bvm ], 

Jk!, Jr? 


and 


.CIIF'lli” 


(5.3) 


(5.4) 


E[ J' J \u E (0,y)-u 0 (x)\g s (x-y)dxdy\ < e[ J' J« e (0, x) - u 0 (x)\ dx\ + 5£[[m 0 Ibv(R‘')] (5-5) 

We combine (15.4b and (15.5b in (15.3b to conclude 
E \ J Ju e (t,y) - u (t,y)\dy\ 

'{£[ J' \u e ( 0 ,y ) - Mo(y)|c/y] + c(l + £’[|Molfiv(R‘ i )])- + < 5 ^[l M olfiv(K rf )]} 
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+ Ce c ^ F " l “'(l + £[|MolBy(R<')])ef + Ce + 5£[|moIsv(R‘')]- 
We choose 6 - in (15.6b . and conclude that, for a.e t > 0, 

E\ I \u e (t,x) — u(t,x)\dx\ 

1 Jr* 1 

<C(r){e 2 (1 + £ , [|wolsv(R rf )])( 1 + t) + E il |« e (0, x) - Mo(x)j c/x]}, 
for some constant C(T) > 0, independent of £ , [|wol,sv(R' i )]- This completes the proof. 


(5.6) 


6. Fractional BV Estimates 

In this section, we consider a more general class of stochastic balance laws driven by Levy noise of the 
type 

{ du(t, x) + di v x F(u(t, x)) dt = q(x, u(t , x); z) N(dz, dt ), x€ll 7 -, 

u{ 0, x) = uo(x), x G M. d , 

Observe that, the noise coefficient z/(x, m; z) depends explicitly on the spatial position x. Moreover, we 
assume that ;/(x, u;z) satisfies the following assumptions: 

(B.l) There exist positive constants K > 0 and A* G (0,1) such that 

\tj(x, w;z) - f](y,v,z)\ < ( A*\u - v| + A"|x - y|)(|z| A 1), for all u, v G R; z G R; x,y G R"”. 

(B.2) There exists a non-negative function g(x) G L°°(R d ) n L 2 (R rf ) such that 

|/ 7 (x, m;z) | < g(x)(l + |m|)(|z[ a 1), for all (x, u,z) g R rf x R x R. 


Clearly, our continuous dependence estimate is not applicable for problems of type (16.1b . and primary reason 
for that lies in the nonavailability of BV estimate here. We refer to [(3, Section 2] for discussion on this point 
for diffusion driven balance laws. However, it is possible to obtain a fractional BV estimate. To that context, 
drawing primary motivation from the discussions in @], we intend to show that a uniform fractional BV 
estimate can be obtained for the solution of the regularized stochastic parabolic problem given by 


du £ (t, x) + di v x F e (u e (t, x)) dt = 


I r/ £ (x, u e (t, x); z)N(dz, dt) + eA xx u e (t , x) dt, 
J|zl>o 


( 6 . 2 ) 


where F e , q, satisfy (12.4b . Regarding equation (16.21) . we mention that existence and regularity of the solution 
to the problem (16.21) has been studied in 0]. We start with a deterministic lemma, related to the estimation 
of the modulus of continuity of a given integrable function, and also an useful link between Sobolev and 
Besov spaces. In fact, we have the following lemma, a proof of which can be found in JjJ Lemma 2], 


Lemma 6.1. Let h : M. d —> R be a given integrable function, 0 < <p G C“(R £/ ) and be a sequence 

of symmetric modifiers, i.e., Js(x) = Jj7(y), 0 < / G C“(R), supp(J) c [—1,1], J(~) = /(•) and f J — 1. 
Then 


(a) For r, s G (0,1) with r < s, there exists a finite constant C\ — C\(J, d, r, s) such that 



\h(x + z) - h(x - z)\Js{z)f(.x) dxdz 


<C i 6 r sup |z|" 
|z|<<5 


f 


| h(x + z) - h(x - z)\4>{x) dx. 


(6.3) 
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(b) For r, s e (0,1) with r < s, there exists a finite constant C 2 — C 2 O, d, r, s) such that 
sup I \h(x + z) - h(x)\c/)(x) dx 

|z|<c5 Jr? 

< C 2 d r sup 6~ s f f \h(x + z)-h(x-z)\Js(z)<p(x)dxdz+C 2 d r \\h\\ L i ( wi ) . (6.4) 

o<5<i Jr? Js.i 

Now we are in a position to state and prove a theorem regarding fractional BV estimation of solutions of 

(IQ . 

Theorem 6.2 (Fractional B V estimate). Let the assumptions |(A.1)| |(A.2)| |(B. 1 )| |(B.2)| and \(AA)\ hold. Let 
u e be a solution of Ol with the initial data uq(x) belongs to the Besov space ^( IKJ ) for some p e (^, 1). 
Moreover, we assume that F" € L°°. Then, for fixed T > 0 and R > 0, there exits a constant C(T,R), 
independent of e, such that for any 0 < t < T, 

supftj I \u e (t, x + y) - u e (t, x)\dx\ < C(T,R) 6 r , 

|j|<<5 JxeKR 

for some r 6 (0, j) and Kr := {x : \x\ < R}. 

Proof. Let 0 < cp(x) e C~(M. d ) be any test function such that |V0(x)| < Ctf. >(x) and |A0(x)| < Ccf>{x ) for some 
constant C > 0. Let J s be a sequence of modifier in as mentioned in Lemma 16.11 Consider the test 
function 

Sutracting two solutions uf t, x), u e {t,y) of (16.2b . and applying Ito-Levy formula to that resulting equations, 
we obtain 


fe(u e (t, x) - u e (t,y)) -fe(u e ( 0, x) - u e (0,y)) 

= J' J3f(u e (s,x) - u e (s,y))(div y F e (u e (s,y)) - d\v x Ffufs,x))'jds 

+e x) ~ x) ■ A yyU e (■ r,y))dr 

+ (1 -P)flf(u e (r,x) - u e (r,y) + p(rj e (x, u £ (r, x); z) - rj e (y, u e (r, y); z))) 

Jr= 0 J|:|>0 Jp=0 

X \p e (x, u e (r, x); z) - T] f (y, u e (r, y); z)|“ dp v(dz)dr 

+ \MuAt X) - u e {r, y) + ij f (x, u e (r, x);z) - rj e (y, ufr, y); z)) 

Jr=0 J|z|>0 

- f^(u E (r,x) - u e (r, y))] N(dz, dr). 


To this end, we see that 

PfU — V)(A XX U — AyyV) — (A,, + 2V* ■ Vy + Ayy'jP^U ~ V) ~~ [f (it ~ V)| V X U — VyV| . (6.5) 

Moreover, a simple calculation reveals that 

( a ** + 2V. r • Vy + Ayy)^ tf (X,y) = J S (^f) , 

(v A + Vy )Mx,y) = V0(^)J 5 (^). 
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Using convexity of /3 f and (16.5b . we have 

I I fS({u € (t, x) - u f (t, y))i// s (x, y ) dx dy - I I fe(u € (0,x) - u € (0, y))if/g(x, y) dx dy 

Jr; Jr; Jr; Jwi 


r =0 JR" JR" 


a x + y x — y 

F r 6 (u £ (s,x),u £ (s,y)) ■ V <p{ —— )Js( —— )dxdyds 
' 2 2 

F £ (u £ (s,y),u £ (s,x)) - F £ (u £ (s,x),u £ (s,y)fj ■ V v i//s(x,y) dxdy ds 

x — y x y 

ePf(u € (r,x) - u e (r,y))J 6 {—^-)^(—j-)dxdydr 


)r =0 J\z\>0 


x J P =o 


{u e (r, x) - u e (r, y) + p(r/ e (x , u € (r , x);z) - n e (y , u e (r, y); z))) 


)r =0 J|z|>0 


Notice that since \F f 


x |t 7 £ (x, M e (r, x);z) - r/ f (y, u £ {r,y)\ z)|Va(x,;y) dpdxdy v(dz) dr 
(u f (r, x) - u e (r, y) + rj £ (x, u £ (r, x); z) - rj £ (y, u e {r, y); z)) 

- Pt(u £ (r, x) - u £ (r,y))]ifrg(x,y) dx dy N(dz, dr). 

?(m, v) - F £ (v, u )| < CIIF'Hk, %\u - v|, we obtain 

e\ I I pAu e (t,x) - u £ (t,y))i//g(x,y)dxdy] - E\ I I pAu e (0,x) - u e (0,y))i/fs(x,y)dxdy\ 
l jr;jr; J l jr;Jr; 1 

< CIIF'IU f e\ f f \u £ (s,x) - u £ (s,y)\(f>(^-^)J s (^——)dxdy\ds 
J s-0 l Jr;Jr; 1 2 2 J 


+ C\\F'\\^E 
+ Q\F'\\o4E 
+ Ce 


h =o Jr; 


)s=o jr; 


+ E 


U £ (s,x ) - u £ (s,y)\<pC^-^)Js(^-^-)dxdyds] 

u £ (s,x ) - M e (s,y)|0( -^-)|V v / 5 ( : ^p-)l c/xafyds] 

f4f f |w e (r, jc) - u £ (r,y)\j s (^-—^-)<p(^^-)dxdy\dr 
J r = o Jr; Jr; 2 2 j 

(«e(r, x) - u £ (r, y) + p(tj £ (x, u £ (r, x); z) - rj £ (y, u £ (r, y); z))) 

x \r/ £ (x, u £ (r. x);z) - rj £ (y, u e (r,y); z)\~i/fg(x,y) dp dx dy v(dz)dr], (6.6) 


r= 0 J|:|>0 Jr; Jr; Jp= 0 


where we have used |A0(x)| < Ccf >(x). 

As before, one can use Cauchy-Schwartz inequality along with uniform moment estimate (12.5b to con¬ 
clude 


C\\F'\UE 
+ C\\F'\\ aa £E 


u £ (s,x ) - u e (s, y) | (pi ~Y ~) dx dy <f.v] 

u £ (s, x) - KeCs, y)|^( )|V,4( -z^-)\ dx dy c/s] 


<C||F'|| 00 (f+|)||^|| £ - (Ri) 4. 
0 


(6.7) 
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Next, we focus on the last term of (16.61) . To estimate that term, we first let 


a = u e (t, x ) - u e (t, y ) and b = t) £ {x, u £ (t, x); z) - rj £ {y, u e {t, y); z). 


Observe that 


b 2 ft(a +pb) = ( rj £ (x,u £ (t,x);z ) - r] £ (y,u e (t,y);z)) 2 ft (a + pb) 


< (| u £ (t,x) - u e (t,y)\ 2 + K 2 \x — y| 2 )(l A \z\ 2 )ft(a + pb) 
= (a 2 + K 2 \x - y\ 2 )ft{a + p b) (1 A |z| 2 ). 

As before (cf. 13.5b . one can use assumption |(B.l )| on r/(x , n; z) to conclude 

0 < a < (1 — T*) _1 (fl + pb + K\x — y|). 

In view of (16.8b . we have 

K\x - y| 2 

\x-y\ 2 


( 6 . 8 ) 


b l ft(a +pb)<(l- AT (a +pb + K\x- y\) ft (a + pb) (\z\ z A 1) ■ 


(Izr A 1) 


and hence 
E 


)r= 0 J|:|>0 

< E 

< 


)r= 0 J|z|>0 JR^ J Ri 

C^+ 5 -)tU{-)\T 


< 2(1 - AT\a + pbfftia + pb)(\z\ 2 A 1) + C(K, A *)——(|z | 2 A 1) 

< [2(1 - AT 2 C{ + C(K , T*) |A ~- Vl ](\z\ 2 A 1), 


b 2 ft (a + p b)tf/ s {x, y) dp dx dy v(dz) dr J 


■ AT 2 C% + C(K,A*) ^ X £ }(|z [ 2 A 1 )i// s (x,y)dxdyv(dz)dr\ 


i j p =o 


(6.9) 


Now we make use of (12.9b . ( 16.7b to (16.9b in ( 16.6b and conclude 

E I I \u £ (t,x) - u £ (t,y)\jft^—- )4>( ^--)dxdy 
Jr? Jr? 1 2 2 

<E f f |« e (0,x) - Me(0,y)|/5(^r -)<K— 7 T-)dxdy 

. Jr? Jr? i i 

+ C(1 + HF'Hoo) f E\ f f I u e {s,x)~ u e {s,y)W^-^)J s {^—^)dxdy\ds 
J s = o 1 2 2 J 

£ 6 2 

+ CIIT'Hoo^ + -)||</>||l°°(r<() Vr + C(£ + ( 1 ll0(-)lli“(M d ) + 

A simple application of Gronwall’s inequality reveals that 

n \u e {t,x) - u e (t,y)\js(-—r-^-)(p(—^-)dxdy 

, i 11 

exp (t C(1 + in| L ~))£[ \u e (0,x) - uMy)\j S (^ l m^ l )dxd } ] 
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( 6 . 10 ) 


+ exp (f c(i + iniL»))[c(ini 00 (^+ |)Vi+(f + j)t)\m L ~m +c t\w\LW)\- 


Chosing f = C5? in (16. 1 Ob . we obtain 

E \ f f |w e (f,Jf) - u £ {t,y)\j s {^—— dx dy] 

L Jr? Jr? 1 2 2 J 

< C(r)£[ f f |m<t(0,x) - M e (0,y)|y a ( : ^— )0(^-~)dxdy\ 

L Jr* Jr* 1 2 2 J 

+ C(T){(6- + Vb)||0|| L°°(R*) + 5 2 II0IIl 1 (R-'))> 

for some constant C(T ) > 0, independent of e. 

Now we make use of the following change of variables 

x — y , x + y 
* = —' and ^— • 

to rewrite the above inequlity (dropping the bar). The result is 
*iu \u e (t, x + y) - u e (t, x - y)\Js(y)<p(x) dx dy J 

< C(T)E UJ) u £ ( 0, x + y) - u £ ( 0, x - y)|y^(y)0(x) dx Jy] 

+ C{T)(i6 2 + Vb)||0||/,«.(R</) + b 2 ||0||ii(Rd)) 

In view of (16.4b of the Lemina [6~il we obtain for r < - 

sup£ I \u £ {t,x + y) —u £ {t,x)\<p{x)dx\ 
l3i|<# l Jr? j 

< C 2 b' sup b”5£j I I \u £ {t,x + y) - u £ {t,x - y)\js{y)4>{x)dxdy\ 

o«y<i L Jr? Jr? j 

+ C2<)'E[\\u £ (t, OIIl^r^)]- 

Again, by (16.3b of the Lemma l6Tl and by (16.1 It . we see that for r - ) and s > ^ 

sup <5^2 £■[ I I \u £ (t,x + y) - u £ (t,x - y)\js(y)c/>(x)dxdy] 

o<<s<i l Jr?Jr? j 

< C(!T) sup b~ 2 /?[ I I \u £ (0,x + y) - u £ (0,x - y)\j$(y)(/)(x)dxdy\ 

o<cS<i L J r? Jr? j 

+ C(7’)(||0|| z 

< C(T)Ci supflyr^fif ( \u £ (0,x+ y) - u £ (0,x)\(p(x)dx\ 

|y|<rf \ ’ L Jr* j 

+ C(r)(||^|| z ,«. (K . 

< C(T) e [\\ u 0 \\ b ^ a (R<')] II0 IIl”(R-') + C(7’)(||0IIl“(r<o + ||0||z.i(R rf ))- 

Now we combine (16.12b and (16.13b to obtain 

sup E\ I | u e {t, x + y) - u £ (t, x)|^(x) dx\ 

|y|<5 L JR? J 


( 6 . 11 ) 


( 6 . 12 ) 


( 6 . 13 ) 
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< C(T)6‘ 


^[llWoll^JR^)] + 1 jll0lli“(R'') + II0 IIl 1 (R‘O 


+ C 2 6 r E 


I u c(t, ■)llz, 1 (R rf ) 


Let Kr — {x : \x\ < R}. Choose f 6 C“(R fl ') such that <p(x) — 1 on Kr. Then, for r < i, we have 


sup E 

|v|<(5 


I | u e (t, x + y) — lift, x)| clx 
•Jk r 


< C(T, R)6 r , 


which completes the proof. □ 

In view of the well-posedness results from Q], we can finally claim the existence of entropy solutions 
for (16.11) that satisfies the fractional BV estimate in Theorem 16.21 In other words, we have the following 
theorem. 

Theorem 6.3. Suppose that the assumptions |(A.2)| |(A.3)| |(A.4)| |( B. 1 )| and |(B.2 )| hold and the initial 
data uq belong to the Besov space B^ oo (M. d )for some p e (|, 1) and 




< oo, for p = 1,2, • 


(6.14) 


(a) Then given initial data uq, there exists an entropy solution of (IQt such that for any t > 0, 

E[ \\u(t, ] < oo, forp = 1,2,- •• . 

Moreover, there exists a constant Cj such that, for any 0 < t < T, 


sup£ 

|v|<J 


I \u{t, x + y) - u{t, x)\ dx 
Jk,, 


< Cj 5 r , 


for some r e (0, 4) and Kr := {x : \x\ < R}. 

(b) Let the initial data uq only satisfies (16.141) . Then there exists an entropy solution of dot such that for 
any t > 0, 


£[ II u(t, •)ll^ (Rrf) ] < oo, for p = 1,2, • • • . 
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